Diffuse inverse Compton and synchrotron emission 
from dark matter annihilations in galactic satellites 
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Annihilating dark matter particles produce roughly as much power in electrons and positrons as 
in gamma ray photons. The charged particles lose essentially all of their energy to inverse Compton 
and synchrotron processes in the galactic environment. We discuss the diffuse signature of dark 
matter annihilations in satellites of the Milky Way (which may be optically dark with few or no 
stars), providing a tail of emission trailing the satellite in its orbit. Inverse Compton processes 
provide X-rays and gamma rays, and synchrotron emission at radio wavelengths might be seen. We 
discuss the possibility of detecting these signals with current and future observations, in particular 
EGRET and GLAST for the gamma rays. 
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I. INTRODUCTION 



> , 

It is almost universally accepted that most of the matter in the universe is non-baryonic. This dark matter is the 
s ! . chief constituent of gravitationally bound objects from dwarf galaxy scales and larger. Identifying the nature of dark 
' matter is one of the most important problems in astrophysics, cosmology, and particle physics. 

Perhaps the best motivated candidate for cold dark matter is the lightest of the so-called neutralinos arising in 
supersymmetric extensions to the standard model p|. These are the spin- 1/2 Majorana fermion counterparts of the 
neutral gauge and Higgs bosons, and are expected to have masses at the weak scale (of order 100 GeV). This scale is 
intriguing as the relic density of a stable particle with this mass and corresponding cross section turns out to be of 
order the critical density, as we observe the matter density to be today. Any stable particles with weak scale masses 
could thus naturally account for the dark matter. In this paper, we will focus on supersymmetry, but our conclusions 



are fairly generic to dark matter candidates at the TeV scale. 



We outline a new signature of annihilating dark matter in satellites of the Milky Way galaxy. Prospects for 



detecting high energy photons as dark matter annihilation products, primarily from the ir° decays that are generic 



to hadronization processes, have been discussed for many years (for a sample see Ref. 0). Necessarily coming with 
these photons are high energy electrons and positrons from the analogous tt^ decay chains. Charged particles suffer 
complicated motions in the galactic magnetic field, and furthermore they lose energy to synchrotron and inverse 
Compton processes. Searching for the synchrotron emission from the galactic center 0] and from galactic satellites 
has been discussed previously, though neglecting the diffusion of the charged particles. We will show that the inverse 
Compton emission, extended over a large area from the charged particle annihilation products may be observable for 
some models of particle dark matter and of galactic satellites. 

II. SUPERSYMMETRIC CANDIDATES 
A. Particle physics model 

In the Minimal Supersymmetric Standard Model (MSSM) the lightest of the superpartners (LSP) is often the 
lightest neutralino. The latter is a superposition of the superpartners of the neutral gauge and Higgs bosons, 

Xx = N U B + Ni 2 W 3 +N 13 H^+N U Hl (1) 
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FIG. 1: Spectra of photons and leptons from annihilations in two characteristic models. The solid curves depict the combined 
electron and positron spectra, and the dotted curves show the photon spectra. Both of these models have acceptable values 
of ft x h 2 . The model on the left has an appreciable Higgsino fraction, and thus annihilates efficiently into W pairs. The direct 
decays W — *■ e ± u are responsible for the shelf in the lepton spectrum at energies right below maximum. The model on the 
right is a fairly pure gaugino, thus the final state is nearly always qq, and the spectrum is featureless. 



With R-parity conserved, this lightest superpartner is stable. For significant regions of the MSSM parameter space, 
the relic density of the stable neutralino is of the order £l x h 2 ~ 0.1, thus constituting an important (and perhaps 
exclusive) part of the cold dark matter. Note that Q x is the neutralino density in units of the critical density and h 
is the present Hubble constant in units of 100 km s _1 Mpc -1 . Current observations, included those of the WMAP 
satellite and the Sloan Digital Sky Survey (we take the WMAP values), favor h = 0.7lig gg and a matter density 
tt M h 2 = 0.135lg;oog, of which baryons contribute a small amount tt B h 2 = 0.0224 ± 0.0009 [1 @ If we assume that 
neutralinos are the only constituent of matter (neutrinos are a small contribution, fl u h 2 < 0.0076), we can then infer 
£l x h 2 — 0.113lgQQg. We will apply a generous 3a constraint on the relic abundance from the WMAP data as follows: 
0.086 < n x h 2 < 0.137. 

Using the DarkSUSY code l7|, we have explored the supersymmetric parameter space in both a phenomenological 
MSSM H H [H [H E E El El and in minimal supergravity (using the IS A JET code El). Each model is 
subjected to current accelerator constraints on masses of superpartners and Higgs bosons 0, Il8| and on the b — > sj 
branching ratio fl9l . Crucial for studies of dark matter, the relic abundance of neutralinos Q x h 2 is calculated based 
onRefs. BE IDT 



B. Photon and electron spectra of annihilations 

The only detailed information about models we require for this study (apart from the relic density constraint) is the 
annihilation cross section at non-relativistic velocities (av) (accounting for the factor of two due to identical particles 
in the initial state) and the spectrum of annihilation products, in particular the photons, electrons and positrons. 
The photons are primarily from the 7r° decays associated with hadronic final states, whereas the leptons can be from 
the associated decays, or more directly from massive gauge boson decays. In Fig. ^ we illustrate two models in 
the generic MSSM, one with a typical featureless spectrum of annihilation products, and the other clearly exhibiting 
the feature from H /± decays. 
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III. CHARGED PARTICLE PROPAGATION AND DIFFUSE EMISSION 



The calculation of the diffuse emission from the annihilations proceeds in several steps. First, the time dependent 
density of charged particles is determined according to a diffusion model. The charged particles are trapped by the 
galactic magnetic field, which extends for several kpc from the stellar disk. Here we stress that we are concerned 
only with those galactic satellites that are currently (or in the recent past) within this "diffusion zone". A time 
dependent treatment is necessary because a satellite moving with the typical galactic velocity of 300 km s _1 crosses 
the diffusion zone in roughly one diffusion time: both timescales are of order tens of millions of years. The second 
step in the calculation is to calculate the column density of particles along various lines of sight, as a function of 
particle energy. Lastly, the inverse Compton and synchrotron spectra can be calculated from the particle spectrum 
under some assumptions about the galactic magnetic field and radiation fields. 



A. Diffusion model 



The propagation of charged particles in the tangled galactic magnetic field can be modeled as diffusion. Electrons 
and positrons lose energy rapidly to synchrotron radiation, and also to inverse Compton scattering on the cosmic 
microwave background (CMB) and on starlight. We will use the diffusion model of Ref. ^ij. More sophisticated 
models are available, e.g. the semianalytic treatment of Ref. [2]| and the fully numerical GALPROP model , but 
our model is simple to implement, yielding quantitatively similar results, and affords an intuitive understanding. We 
will make the assumption that the charged particle velocity distribution is locally isotropic, namely that the particles 
"stick" to the galactic magnetic field, and that the velocity distribution contains no record of the initial conditions. 
This is important as we consider sources (satellite halos) that move a significant distance in a diffusion time. If the 
velocity distribution was isotropic in the satellite frame (as could be the case if the satellite had a significant magnetic 
field of its own), our results would be different. 

The time dependent diffusion-loss equation in homogeneous space for the density of charged particles as a function 
of energy dn/dE is the following, 

»* 1^*+ 1 » (2) 



dtdE v ' dE E t dE V dE 

where K(E) = K Q 3 a [l + (22/ 3 GeV)"] is the energy dependent diffusion constant, with K = 3 x 10 27 cm 2 s _1 and 
a — 0.6. Energy loss due to any electromagnetic process where the exchanged energy is much less than m e c 2 scales 
as E 2 , and is proportional to the energy density in photons (or equivalently magnetic field). We consider synchrotron 
losses dues to a 3 /iG galactic magnetic field (0.2 eV cm" 3 ), inverse Compton losses from the CMB (0.3 eV cm -3 ), 
and inverse Compton losses due to starlight (0.6 eV cm -3 ) |2g. We thus set dE/dt = E 2 /(E t), with E = 1 GeV, 
and t = 10 16 s. Lastly, Q is the source function. Solving instead for F — E 2 dn/dE, we find 

f d K(E) , 1 d \ 

EPdt E 2 ^ ~ E^t'Be) ^ = < 3 > 

We can calculate the Green function for this operator simply by 4-dimensional Fourier transform in x — > k and t lu, 

CG = 8 z (x-x')5(t~t')5{E-E') -» (4) 



_L Hu + R - (£)tV _i__| 



G = -l—e^'+^'HiE-E'). (5) 
(2tt) 4 y ' w 



We should note that the Green function vanishes for t' > t and for E' < E because time increases monotonically, and 
energy decreases monotonically. This equation is easily solved for E ^ E'\ 

G = Q oe i l UJT ( u - u ')- K o Tk2 ( v ~ v ')] f (g) 

with u — Eq/E and v — 3 a u+u( 1 ~° 1 ' /{I — a). Hereafter we will denote the difference X — X 1 = JS.X for some variable 
X. We now apply a jump condition to fix Go: for E > E' . G = 0. The condition is 

-^=^^=(2^ eiM ' + ^' W 
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and we thus derive the Fourier transform of the Green function, 

fi _ ^0 T i(u>(TAu+t')+k-x'-K Tk 2 Av) (o\ 

~ (2tt) 4 ' [ ' 

Inverting the transform, we find the expected behavior that time and energy propagate in lockstep, yielding the 
free-space Green function 

Gfree = tT**l* S(At - T Au), (9) 

defining the effective diffusion length D 2 = 4KqtAv — (3.550 kpc) 2 At>. This scale indicates that particles travel only 
a few kpc before losing most of their energy (Av ~ 1). 

The galactic magnetic field has a limited extent, which can be modeled simply by defining a diffusion zone at the 
boundary of which particles freely escape. Diffusion models typically require that the height of the diffusion zone is 
larger than L = 3 kpc from the disk on both sides, and the radius is at least 20 kpc, and probably larger. For example 
the best GALPROP models |22j use L = 4 kpc and a radius of 30 kpc. As the radial boundary is far from the Earth, 
we can safely neglect it in considering satellites within 10 kpc of us. We thus model the diffusion zone as an infinite 
slab of thickness 2L, and we take L = 3 kpc. We impose the boundary conditions that the density vanish at z = ±L, 
which can be effected by a series of image charges at positions x n = x, y n = y, z n = (— l) n z + 2Ln, 

oo 

G 2L (x,x') = J2 (-l) n <W2,0 (10) 

71— ■ — OO 

The density of charged particles is now derived, 

rl<r\ P POO 

F = E 2 —= d 2 x' dz' df dE' G 2 L{x,x',t,t',E,E')Q(x',t',E'). (11) 
dE J R 2 J_ L J_ 00 J E 

B. Diffuse emission: total power 

For diffuse emission we are interested not in the local density of particles, but in the column depth of particles, 

E 2 ^=\dlE 2 ^. (12) 
dE J dE V ' 

Placing the observer at the origin, at an angle 8 from the galactic plane and an angle a from the line y = 0, we find 
that for a distance I from the observer, 

x = £ cos 6 cos a, (13) 
y = £ cos S sin a, (14) 
z = £smd. (15) 

For simplicity we define x — xjl. The integral in £ can be applied directly to the free Green function, yielding 
the column depth at coordinates (a, 6) due to a source at x' (and we truncate at the edge of the diffusion zone, 
W = L/\ sin<5|), 



. / ^max X • X . 

erf — — erf 



D J V D 

The Green function satisfying the boundary conditions is now simply 



5{At-rAu). (16) 



G% L (a,5,x')= J2 (-l) n Gl cc (a,S,x' n ), (17) 



and the column density of charged particles is 



E 2 %= f d 2 x' I dz' I dt' [ dE'G 2 T L (a,6,x*,t,t',E,E')Q{$',t',E'). 
dE J R 2 J_ L J_ 00 J E 



(18) 
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We consider a toy model for dark matter clump undergoing annihilations. Assuming a very cuspy profile, most 
annihilations will occur very close to the center, and we can assume the clump is a point source. Taking the clump 
to be at position X = (X, Y, Z) today (t = 0), moving with constant velocity V, we can write the source function as 



Q(x,t,E)=T^5 3 (x-X-Vt), 



(19) 



where T is the annihilation rate in the clump, and d<p/dE is the spectrum of per annihilation. The delta functions 
simplify matters greatly, leaving only an integral in energy. The charged particle density (both e + and e~) at the 
Earth and the column depth are given by 



E 2 



, dn 
dE 



deb EqTt 
dE' (ttD 2 ) 3 / 2 



(20) 



71— — OG 



E^ = 
dE 



dE' 



E Tr 



OG 



dE' 2ttD 2 ^ 

TI— — 00 



(-1)' 



\]/D 2 



erf 



D 



erf 



-x ■ w n 



>(L-\w,\). (21) 



The effective position w = X — VtAu, and w z ,n = (—l) n (Z — V z tAu) + 2Ln. We recover the steady-state solution 
by taking V = 0, indicating a source that has been at location X for all time. 

We observe that looking at lower energies implies looking back in time, to when the source was in a different 
position. This leads to the formation of a "wake" of diffuse emission extending from the current position of the source 
(which will be a bright, concentrated source of gamma rays) in the direction it came from. The spatial signature is not 
unlike a comet, with a bright "nucleus" of gamma rays from the 7r° decays, and an extended tail of diffuse emission. 

The energy loss time r represents all energy loss mechanisms: inverse Compton scattering from both the CMB 
and starlight, and synchrotron emission. The total power emitted in any one of these is inversely proportional to the 
timescale for the individual process: w 2t, tomb ~ 4r, and r sync rs At. Integrating in energy, 



Px = 



1 



dEE^ 



4irE Q Tx J dE' ^ 

and we note this has units power per area per solid angle. In Fig. [21 we illustrate contours of the total power in diffuse 
emission for several geometries of satellites. Finally, we note that the density of charged particles at the position of 
the Earth coming from the clump is typically unobservably small. 



C. Diffuse emission: inverse Compton spectrum 



For electrons of Lorentz factor 7 = E/m e scattering from a blackbody spectrum of photons with hvT = kT. 
photon spectrum is given by the usual inverse Compton formula |23| , integrated over the blackbody spectrum 

" 4 dy2ln(y/A) + l + 4/y-y/2 



1 dP 



x_ 

Px dx 128C(3)~ Jo V exp(x/y) 

with x = is/(vt"/ 2 )- It is easy to show that this expression integrates to unity. This formula is inherently non- 
relativistic, implying that jhi'T <C m e c 2 . For the CMB, this is easily satisfied, as hvx = 0.2348 meV for a 2.725 K 
blackbody spectrum. For starlight, we assume a blackbody spectrum with T = 3800 K [24[, thus hvx = 0.33 eV. The 
non-relativistic condition is satisfied for electron energies E < TeV. This is true for most if not all of the annihilation 



9 



1 



the 



(23) 



products we will consider, 
frequency) , 



We use Px = E /(Eqtx) to find the spectrum from a single particle (in power per 



dPx 



dv 128<(3) E t x vt 



dy2\n{y/A) + l + A/y-y/2 



y 



exp(x/y) - 1 



The spectrum from the full column of particles is then simply (in power per area per frequency) 



dP x = m 2 e 
dv AkEqtxVt 



dE 



da 
dE 



1 dP x 
Px dx 



E 2 



(24) 



(25) 



Detecting a diffuse signal in gamma rays will necessarily be hampered by at least the extragalactic background. 
The EGRET satellite measured this to be El 



E 



dE 



(3.30 ±0.15) x 10" 



E 



0.451 GeV 



-1.10±0.03 



sr 



(26) 
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longitude (y) (degrees) longitude (y) (degrees) 




FIG. 2: Total power in diffuse emission. The normalization is arbitrary; this figure is meant to illustrate the spatial extent of 
the signal. The contours are separated by 0.25 decade in power. Four satellite geometries are illustrated, each with a velocity of 
300 km s _1 . The annihilation spectrum is that of the left panel of Fig. Q In each case the current angular position is indicated 
with a cross. Top left: (x, y, z) = (8, 0, 1) kpc, v = (0, 5/13, 12/13). In this case the clump has crossed more than half of the 
diffusion zone. Top right: (x,y,z) — (12,0,4) kpc, v = (1, — l,l)/v3- This illustrates the case where the clump has passed 
completely through the diffusion zone, leaving behind a wake that will dissipate shortly. A hint of the edge of the diffusion zone 
can be discerned at latitude ~ 14°. Bottom left: (x,y,z) = (8,0,2) kpc, v = (0, 1, — 1)/\/2. In this case the clump has recently 
entered the diffusion zone, thus the wake has not had much chance to form, and the diffuse flux drops very rapidly. Bottom 
right: (x, y, z) — (8, 2, 0) kpc, v = (0, 12, — 5)/13. This final panel depicts a clump moving nearly parallel to the galactic plane, 
thus the time spent inside the diffusion zone is relatively long, and the tail becomes more extended. 
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D. Diffuse emission: synchrotron spectrum 



The total power in synchrotron radiation is similar to that in inverse Compton emission in that both are proportional 
to E 2 . The spectra are different however. 

i dp sync 27Vs r f°° . 

p j=^~ x dasma dyK 5/3 (y), (27) 

-fsync ax JO Jx/sina 

and now x = v/(vbJ 2 ), with vb = 3eB / (4nm e c) — 12.6 (.B/3/iG) Hz. This is just the usual synchrotron formula, 
isotropically integrated over the angle the particle trajectory makes with the magnetic field. It is again easy to show 
that this expression integrates to unity. The spectrum from the full column of particles is then easily recovered from 
Ea. 1251 replacing the {dP/dx)/P with the synchrotron formula, and replacing i>t with vb (and we take B = 3/iG). 
The chief diffuse background for this synchrotron signal is the CMB, which peaks at 160 GHz. At 1 GHz, the power 
is 2.09 x 10 4 Jy sr _1 . In Fig. |31 we plot the spectra of inverse Compton and synchrotron emission from a line of 
sight that is 5 degrees behind the current position of a clump. We will neglect complications such as synchrotron 
self-absorption and synchrotron self-Compton processes as we calculate the photon density to be much smaller than 
in previous treatments |4| due to the diffusion of the source electrons. 



IV. GALACTIC HALO SUBSTRUCTURE 



Many authors have discussed the distribution of substructure in galactic halos. There are two important issues 
relating to the detectability of annihilations in these structures, both their mass distribution and their density profile. 
Again we stress that we are interested only in those satellites within the diffusion zone. As the diffusion zone is fairly 
extended, there should be a significant number of such satellites: even out to a radius of 30 kpc, the diffusion zone is 
roughly 15% of the total volume of the galactic halo. 



A. Density profiles 

We first discuss the density profiles of the satellite halos that we hope to detect. Numerical N-body simulations 
find that structures at all scales diverge as a power law at small radii. In particular, the popular NFW profile [26| is 
given by 



p(r) = po (-) -7- 
\r / (a 



(a + rf 



(28) 



These profiles can be normalized according to simulations, based on the virial mass and radius, where the density is 
roughly 340 times the cosmological density (for the concordance model) p c = iVtH^/^TiG w 1.4 x 10 -6 GeV cm -3 . 
The relevant parameter is the concentration C — r v [ r /a, with r v j r = 260(M vir /10 12 Mq) 1 / 3 kpc [2]]], finding 



C 

Po 
M vir 



107 



-0.084 



340C(1 + C*) 2 p c , 



47r Po ( -^r 



ln (1 + c) 



C 



l + C 



The specific annihilation rate is then given by 

T _(av) Po l-(l + C)~ 3 



{av)p 



1 



3m 2 ln(l + C) - C/(l + C) 



3to 2 In C — 1 



(C » 1). 



(29) 
(30) 
(31) 

(32) 



The Moore profile has a steeper central cusp with a divergent total flux |2q . Thus a minimum radius r m j n needs to 
be defined to regularize. The profile is given by 



p{r) 

r 

M vir 



Po 



1.5 



1 



r) [1 + (r/a) 1 - 5 ] ' 



(<Tv)p ln(a/r mill ) 
m 2 In C 



(33) 
(34) 
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FIG. 3: Spectra of diffuse emission in inverse Compton and synchrotron radiation. The clump is located at (a;, y, z) = (8, 0, 1) 
kpc, with a velocity of 300 km s" 1 in the v — (0,5/13, 12/13) direction, as in the top left panel of Fig. The annihilation 
spectrum is that of the left panel of Fig.0 The line of sight illustrated trails the current location of the clump by 5 degrees in the 
—v direction. The total annihilation rate is 10 39 s _1 . Top left: Spectrum of inverse Compton emission from interactions with 
starlight, assuming an energy density of 0.6 eV cm -3 with a 3800K blackbody spectrum. The dotted curve is the diffuse gamma 
ray background. Top right: Inverse Compton on the CMB. Bottom left: Synchrotron emission assuming a 3 fiG magnetic field. 
The dotted curve is the CMB. Bottom right: Column depth of particles as a function of energy, weighted by the energy loss 
rate. It is clear that most of the diffuse power is radiated by relatively high energy particles. 



For cusps p oc r a steeper than a = 1.5, we find the behavior (with the proper limit as a — > 1.5) 



(<Tv)p Q (a/r n 



\2a~3 



M v 



{2a -3) In C 



(35) 
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B. Mass distribution of subhalos 



Simulations indicate that the mass distribution of subhalos is a power law, given by 

with a normalization that about 500 clumps with masses larger than 10 s M Q are found in a halo like that of the Milky 
Way 29]. This means that there will be only a handful of clumps more massive than 10 10 Mq. At M v - U = 1O 9 M0, 
the concentration is C ~ 19. The specific annihilation rate for a clump is easily determined, using (av) p c /m 2 = 
5 x 10 24 ({av) /10~ 27 cm 3 s _1 ) (m/GeV)~ 2 Mq s -1 . Using po ~ 340/9 c C 3 for any of these halo models, we can find a 
simple expression for the total annihilation rate in a "large" clump, defined as M v ; r = 10 9 Af Q , r v ; r = 26 kpc (and 
a = 1.4 kpc): 

r NFW = 2xl0 39 ^s- 1 , (37) 

m GcV 

T a = 4 x 10- (a/7 T )2a : 3 - 1 W*L s~\ (38) 
2a — 3 m GcV 

In supersymmetric models, (av) 27 / 'm^y < 10~ 2 . The Moore profile has a logarithmically diverging flux, but the en- 
hancement factor is less than 40 if annihilations empty the central cusp, and much less than that for other mechanisms. 
For example, if interactions with a central black hole sweep out the central region (of order 1 pc), the enhancement 
is less than 10. This may be the case for the full galactic halo, but for dark clumps without stars, there may be no 
central black hole. Steeper profiles can in principle have much larger enhancement factors, though with a cutoff of 
1 pc, the enhancement is less than of order 1000 even for a = 2. In the next section, we will see that annihilation 
rates of order 10 38 s _1 are needed for point sources to be confidently seen in diffuse emission, thus a > 1.5 is probably 
required. Coincidentally, steeper profiles look more like point sources, so clumps with high enough annihilation rates 
will be easy to model from this perspective. However, the NFW profile is not too far from detectability for a 10 9 M Q 
clump, if the annihilation cross section is the maximum allowed. Some clumping of the halo might tip the balance 
and make even this profile a candidate for detection in inverse Compton emission. Here we note that the effect that 
dumpiness in galactic halos has on dark matter detection has been discussed at length, e.g. Ref . |30| . 



V. UNIDENTIFIED EGRET SOURCES 



The EGRET gamma ray survey discovered a large number of point sources, a large fraction of which remain 
unidentified [3]]]. These sources typically have fluxes of order 

We investigate the possibility that some of these might be Milky Way satellites undergoing annihilations at their 
centers. The annihilation scenario requires that there is no other concentrated emission than the gamma rays. Other 
astrophysical sources of just about any type are expected to have emission at other wavelengths as well. These sources 
(identified or not) are not measured below about 100 MeV, thus the turnover expected for tt° photons at m^o /2 = 67.5 
MeV could not have been observed. However, the spectrum of an annihilation source should be closer to flat around 
100 MeV, as seen from Fig.Q] thus the EGRET point sources are typically not very good candidates for annihilation 
radiation. However, they do set the scale of possible annihilation sources quite well. 

The yield of gamma rays predominantly from ir° decays in the hadronization of annihilation products is typically 
dN/d In E ~ 2 annihilation -1 at energies of 0.1 GeV. This means the inferred rate of annihilations for these sources 
is 5 x 10~ 8 cm -2 s _1 . For a clump at a distance of 8 kpc, the inferred annihilation rate is 4 x 10 38 s" 1 . As we have 
seen in the previous section, this rate likely requires a central density cusp at least as steep as r~ 15 . 

We ask the following question: can GLAST confirm or rule out the hypothesis that some of the unidentified EGRET 
point sources are the annihilating cores of very steep profile dark matter subhalos? We argue that the answer is yes, 
for two reasons. First, GLAST will have improved performance an energies below 100 MeV, down to a threshold of 20 
MeV. Certainly, the turnover in spectrum at about 70 MeV that should be present in an annihilation source should 
be visible. Second, the diffuse inverse Compton emission associated with a few times 10 38 annihilations per second, 
covering hundreds of square degrees, should yield both a significant fraction of the background and a statistically 
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Energy 
(GeV) 


Size (68%) 
(sr) 


S / BG 


BG (EGRET) 
(counts) 


BG (GLAST) 
(counts) 


S/N (EGRET) 


S/N (GLAST) 


0.03 - 0.1 


0.078 


0.551 


250. 


5,630 


8.72 


41.4 


0.1 - 0.3 


0.084 


0.961 


75.2 


4,230 


8.33 


62.5 


0.3 - 1 


0.093 


1.447 


23.6 


2,650 


7.03 


74.6 


1 - 3 


0.102 


2.098 


7.25 


1,630 


5.65 


84.8 


3-10 


0.102 


3.209 


2.06 


463. 


4.60 


69.0 


10 - 30 


0.076 


5.322 


0.43 


97.1 


3.49 


52.4 


30 - 100 


0.042 


8.019 


0.07 


15.1 


2.08 


31.2 


100 - 300 


0.021 


6.562 


0.01 


2.10 


0.63 


9.50 



TABLE I: Detectability of diffuse inverse Compton emission in EGRET and GLAST. An annihilation source of 4 x 10 s 
is assumed, with the spectrum of the left panel of Fig. Q The EGRET exposure is assumed to be 10 9 cm 2 s. For GLAST, we 
assume a 5 year mission, with a total exposure of 2.25 x 10 11 cm 2 s for energies above 1 GeV. In the energy bins below 1 GeV, 
we assume exposures of 50% (0.3-1 GeV), 25% (0.1-0.3 GeV), and 10% (0.03-0.1 GeV) of this value pg]. The size column 
indicates the 68% containment region. The signal / background column represents the total number of photons, inside the 68% 
containment region, compared with the extragalactic background. For a convincing detection, we require that the signal be 
an appreciable fraction of the background, as the anisotropies in the background are in principle considerable. We of course 
also require a significant signal to noise ratio. Lastly, we require a significant total number of photons. Based on these rough 
criteria, the detectability of the diffuse emission by EGRET seems marginal, though GLAST should be able to detect it easily. 
If the annihilation rate were ten times smaller, GLAST should be able to find the signal as well. Note that "3" in the energy 
bands means \/l0, i.e. the energy bands cover 0.5 decade each. 

significant detection. GLAST will have an enormously improved exposure over EGRET: at 1 GeV and above, the 
expectation is 4.5 x 10 10 cm 2 s to any point on the sky, per year j^. In Fig.0] we illustrate the integrated flux from 
inverse Compton on starlight in four energy intervals due to an annihilation source of 4 x 10 38 s" 1 with the spectrum 
of the left panel of Fig.nj We find that such a source should be easy for GLAST to detect, in fact a source ten times 
less active, with 4 x 10 37 annihilations per second should be detectable. These results are summarized in Table. [H 

VI. DIFFUSE SYNCHROTRON EMISSION 

An annihilation rate of 10 39 s _1 produces a significant amount of diffuse synchrotron emission. Previous authors 
have discussed this in the context of clumps [J, neglecting the effects of diffusion. Their calculation essentially 
assumed a point source of synchrotron radiation (actually 10' x 10' f=a 10~ 5 sr), with the same total power as we find. 
We calculate that the diffusion of the charged particles spreads out the synchrotron signal over roughly 0.1 sr, a huge 
area compared with the previous assumption. For this clump, the synchrotron emission exceeds the CMB only below 
1 GHz. In principle, with spectral information the synchrotron and CMB can be untangled, but as the signal lies near 
the galactic plane, this may not be feasible. 

VII. THE GALACTIC CENTER SOURCE 

The galactic center is known to be a bright source of gamma rays, with a broken power law spectrum measured by 
EGRET to be jH 

E % = ( 4 - 2±0 - 02 ) Xl °"(l^lev) " Cm " S "' ( 4 °) 
a = 0.3 ±0.03 {E < 1.9 GeV), (41) 
a = 2.1 ±0.20 (£>1.9GeV). (42) 

We briefly discuss the possibility that this is an annihilation source. At 100 MeV, this spectrum corresponds to 
5 x 10~ 7 annihilations cm~ 2 s -1 . Assuming the galactic center source is at a distance of 8.5 kpc, the inferred 
annihilation rate is 4 x 10 39 s -1 . This is a factor of ten brighter than the bright clumps discussed previously. In Fig.0 
we illustrate the inverse Compton emission that would be associated with the galactic center source. We find that 
the emission corresponds to a small fraction of the diffuse emission near the galactic center [34]], and furthermore, it 
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FIG. 4: Contours of integrated diffuse flux in four energy bands. The clump is located at (a;, y, z) — (8, 0, 1) kpc, with a velocity 
of 300 km s" 1 in the v = (0, 5/13, 12/13) direction, as in the top left panel of Fig. [5] An annihilation source of 4 x 10 s is 
assumed, with the spectrum of the left panel of Fig.Q Two sets of contours are plotted: the solid ones indicate the fraction of 
diffuse flux contained for 10%, 16%, 33%, 50%, 68%, and 95%, while the dotted ones are in units of photons cm -2 s _1 sr _1 , 
separated by 0.5 decade in flux. The flux is very asymmetric at low energy, becoming nearly round at the highest energies. 



falls with distance from the center much more quickly than the observed halo (for an analysis of the gamma-ray halo 
see Ref. (sEJ )■ If the galactic center region were better understood, it might be possible to uncover such a signature. 
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FIG. 5: Contours of integrated diffuse flux in two energy bands associated with the galactic center source. The source is 
assumed to be fixed at (x,y,z) — (8.5,0,0) kpc. An annihilation source of 4 x 10 39 s _1 is assumed, with the spectrum of the 
left panel of Fig.0 The contours are in units of photons cm -2 s" 1 sr" 1 , separated by 0.25 decade in flux. 



VIII. DISCUSSION AND CONCLUSIONS 



We have calculated the spectrum and spatial extent of diffuse emission from the charged particle products of 
dark matter annihilations. In addition to the synchrotron emission discussed previously, we have studied the inverse 
Compton radiation, primarily on starlight photons. We have focused on galactic satellites that are currently within the 
diffusion zone, namely within a few kpc of the stellar disk. For satellites moving with typical galactic halo velocities 
of 300 km s^ 1 , the crossing time of the diffusion zone is of the same order as the diffusion time, thus an inherently 
time-dependent treatment is required. 

For annihilation sources, e.g. galactic satellites at typical distances of 10 kpc, the diffuse emission in both inverse 
Compton and synchrotron extends over roughly 300 square degrees. We have shown that at least in terms of the 
number of photons, the diffuse inverse Compton emission might be detectable by GLAST, assuming bright enough 
annihilation sources. The spatial extent of the emission makes its detection problematic of course. GLAST will 
certainly detect a significant number of point sources in a region of this size. In a future work we will study in detail 
the feasibility of separating these signals. 

As mentioned previously, these results are fairly generic, and do not depend strongly on the particle physics model. 
As we are concerned with the electrons and positrons, we do not even require that the dark matter have hadronic 
interactions. Leptonically interacting dark matter |3a . l37| would still provide photons and electrons, albeit by different 
processes and with different spectra. Such photon sources would be even harder to reconcile with the EGRET point 
sources, but annihilation sources below the EGRET detection limit may be detectable by GLAST in any case. 
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